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THE GLOBAL ATTRACTOR OF THE 2D BOUSSINESQ EQUATIONS 
WITH FRACTIONAL LAPLACIAN IN SUBCRITICAL CASE 

AIMIN HUANG AND WENRU HUO 


Abstract. We prove global well-posedness of strong solutions and existence of the global 
attractor for the 2D Boussinesq system in a periodic channel with fractional Laplacian in 
subcritical case. The analysis reveals a relation between the Laplacian exponent and the 
regularity of the spaces of velocity and temperature. 


1. Introduction 

This paper studies the existence of global attractor for the solution operator S{t) to 
the two-dimensional (2D) incompressible Boussinesq equations with subcritical dissipation. 
The 2D Boussinesq equations read 

dtu -|- u ■ Vw -|- iy{—A)°‘u = —Vvr -|- 9e2, x € D, t > 0, 

, , V • u = 0, X G D, t > 0, 

dt9 + u ■ V6 + k{—A)^6 = /, X G D, t > 0, 

_u(x,0) = uo{x), 6{x,0) = Oo{x), X G D, t > 0, 

where Q = [ 0 , L]^ is the periodic domain, u > 0 the fluid viscosity, and k > 0 the diffusivity; 
u = u{x,t) = {ui{x,t),U 2 {x,t)) denotes the velocity, tt = 7r(x,t) the pressure, 6 = 6{x,t) a 
scalar function which may for instance represents the temperature variation in the content 
of thermal convection, 62 = ( 0 , 1 ) the unit vector in the vertical direction, and / = /(x) a 
time-independent forcing term. Since in this article we consider 2D Boussinesq equations 
with a subcritical dissipation, we assume that the exponents a and /3 satisfy 

( 1 . 2 ) ;3g(^,1). 

Recently, the 2D Boussinesq equations and their fractional generalizations have attracted 
considerable attention due to their physical applications and mathematical significance. 
When a = /3 = 1, the system (1.1) is then called the standard 2D Boussinesq equations, 
which are widely used to model the geophysical flows such as atmospheric fronts and oceanic 
circulation and also play an important role in the study of Rayleigh-Benard convection 
(c.f. [Ped87]). Besides, the 2D Boussinesq system are the two-dimensional models which 
retain the key vortex-stretching mechanism as the 3D Navier-Stokes/Euler equations for 
axisymmetric swirling flows (c.f. [MB02]). 
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In the mathematical respect, the global well-posedness, global regularity of the standard 
2D Boussinesq system as well as the existence of the global attractor have been widely 
studied, see for example [CD80, FMT87, ES94, CN97, MZ97, Wan05, CLR06, Wan07, 
KTWll]. Recently, there are many works devoted to the study of the 2D Boussinesq 
system with partial viscosity, see for example [HL05, Cha06, HK07, DP09, HK09, HKRll] 
in the whole space and [ZhalO, LPZll, HKZ13, Hual5b] in bounded domains. There 
are also many works which considered the global regularity of 2D Boussinesq system with 
fractional diffusion, see for example [WX12, JMWZ14, XX14, YJW14, SW14]. 

In some realistic applications, the variation of the viscosity and diffusivity with the tem¬ 
perature may not be disregarded (see for example [LB96] and references therein) and there 
are many works on this direction, too, see for example [LB96, LB99, SZ13, LPZ13, Hual5] 
where the existence of weak solutions, global regularity, and existence of global attractor 
have been studied. However, the global regularity for the inviscid 2D Boussinesq system 
where = k = 0 is still an outstanding open problem and the study of fractional Boussinesq 
system may shed some light on the inviscid Boussinesq system. To the best of our knowl¬ 
edge, the existence of the global attractor for the 2D Boussinesq equations with fractional 
dissipation has not been addressed yet, which is the goal of this article. 

This work is motivated by the [CC04, Ju05], where the existence of global attractor of 
the 2D subcritical SQG equations has been proved. The key point in [CC04, Ju05] is that 
they proved a positivity lemma for the fractional Laplacian. Armed with the positivity 
lemma, we can similarly obtain the maximum principle for 0 as in [CC04, Ju05] and then 
follow the standard procedure to show the existence of global attractor. We point out that 
our estimates are much more involved and harder than those in [Ju05] since the velocity 
is given by an evolution equation instead of the Riesz potential of 9 and we have to use a 
variant of uniform Gronwall lemma (see Lemma 2.4) to prove the uniform bound for the 
velocity u. 

In this article, we prove the existence of a global attractor for (1.1) in the space x 
where 

(1.3) Si > 2max{l — a, 1 — /3}, S2 > 1, 
and 

(1.4) 0 < S2 — Si < a + fi. 

Condition(1.3) is natural for si in that the smaller a and j3 are, the less regularity we can 
obtain from the equation (1.1)3, and hence we need to assume more regularity for the initial 
data 00 (that is larger si) to compensate for reduced smoothing effect. We need at least 
H^-regularity (that is S2 > 1) for u to show the L^-estimate for u, see Section 3.3. An 
interpretation of condition (1.4) is that the interplay of 9 and u in (1.1) restricts the gap 
between si and S 2 - 

The roadmap of this article is as follows. In Section 2, we introduce the notation, some 
preliminary results, state our main results, invoke the results from [Tem88] to prove the 
existence of the global attractor in x space for si and S 2 satisfying (1.3) and (1.4). 
Section 3 is devoted to the proof of the H^-uniform bounds and for {9,u), where we first 
show the uniform L^-estimates for p > 2 in subsections 3.1-3.3. Then, in subsection 3.4, 
we prove that the compactness of the absorbing ball in x H^'^, which is the key of 
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proof of the compact global attractor. In Section 4, we show that the solution operator 
{S{t) ; t > 0} is continuous from M+ to x and it is continuous from x to 
X fQj- fixed t > 0. These also give a proof of the existence and the uniqueness of 
the strong solution for (1.1). 


2. Notations and preliminaries 

2.1. Notations and function spaces. We denote by C{I,X) the space of all continuous 
functions from the interval I to some normed space X, and by U’{n) the space of the 
pth-power integrable functions normed by 

\\f\\Li>={j \f{x)\Pdx]'’, ||/||l°o = esssup|/(x)|. 

\Jq, J xgq 

We also denote by LP(0, T; X) the space of all measurable functions u : [0, T] i-A X with the 
norm 

ll^ll^yPfO T'A') / ll^llL°°(0,T;Ar) SUp ||'u||jV’ 

' Jo iG[0,T] 

The Fourier transform / of a tempered distribution f{x) on the periodic domain H is defined 
as 

where k = (fci,/c 2 ) is a tuple consisting two integers. We denote the square root of the 
Laplacian (—A) 2 by A and we have 

Af{k) = \k\f{k), 

where |A:| = -s/kf + . We define the fractional Laplacian A^f for s € M by its Fourier 

series 

A7 := 

ke.7? 

For any tempered distribution / on H and s G M, we define the norm 

ii/iIh. = iiA'/iiis = ( E 

\fcGZ2 

and H^{Q) denotes the space of the functions / such that ||/||_fis is finite. For 1 < p < 00 
and s G M, the space consists of the functions / such that / = A~^g for some 

g G The Lf^’^-norm of / is defined by 



IIA^IIlp. 
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2.2. Some preliminary results. We first recall the Gagliardo-Nirenberg interpolation 
inequality, which is used frequently in this article. 

Lemma 2.1 (The interpolation inequality). For any si < s < S 2 , we have 

IIA'slli. <||A*‘9l|*.||A*>9||^‘. 

where s = (5si + (1 — S)s 2 for some 0 < 5 < 1. 

Next, we recall the Gronwall Lemma and Uniform Gronwall Lemma. 

Lemma 2.2 (Gronwall Lemma). Let g, h, y and ^ be locally integrable functions on 
(to, +oo) such that 

d?/(t) 


dt 


< g{t)yit)+ h{t), vt>to. 


then y(t) satisfies the following inequality, 

yit) < y(to)exp ( / g{s)ds ) + / /i(s)exp ( / c/(r)dT ) ds. 

\Jto J Jto \Js J 

Lemma 2.3 (Uniform Gronwall Lemma). Let g, h and y be non-negative locally integrable 
functions on (to, +oo) such that 

dt/(t) 


dt 


< g{t)y{t) + h{t), Vt > to. 


and 


ft+r 


/ t+r 

y{s)d 


s < 03, 


vt > to, 


Vt > to- 


rt-tr 

/ g{s)ds < oi, / h{s)ds < 02 , 

Jt Jt 

where r, 01,02 and 03 are positive constants. Then 

y{t + r) < + 02 ) e“S 

For the proof of uniform Gronwall Lemma, one can refer to [Tem88, pp. 91] and [FP67] 
for its proof. We now present a variant of uniform Gronwall lemma, which is the key to prove 
the uniform estimates for the velocity u. We refer the interesting readers to [GPZ09, Patll] 
for other variants of the uniform Gronwall lemma. 

Lemma 2.4. Let X > 0 be a positive constant and g{t), y{t) be non-negative locally inte¬ 
grable functions on (0, + 00 ) such that 

dy(t) 


and 


I 


dt 

t+i 


+ Xy{t) < g{t), Vt > to. 


5(s)ds < ai, Vt > to > 0. 

Then, there exists t*(to,yo) > 0 large enough, such that for all t > t*(to,?/o); we have 

y{t) < c, 

for some constant C > 0 independent of the time t and initial data yo- 
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Proof. Multiplying and integrating on both side of the equation above from to to some 
t > to, we have 

' ' is. 


Jtn 


*0 

Suppose to+ rn<t< to+ m + l for some integer m, then 

™ rto+k+l 




<e-Ah-io)y(io) + e-^*^( 




e^^5(s)ds 


k=0 


rto-\-K 

J tQ-\-k 


g{s)ds 


( 2 . 1 ) 


< ^ ^m+k+i) 


/c=0 




A - 1 




,2A 


eA - 1 


□ 


We thus finished the proof of Lemma 2.4. 

We now recall an improved positivity lemma due to [Ju05]. 

Lemma 2.5 (Improved Positivity Lemma). Suppose s G [0,2] and 9, A^6 € Then 

[ \6\P-^9A^edx>- [ fAfjeji'j^dx. 

Jn P Jn^ A 

We will use the following Kate-Ponce and commutator inequalities from [KP88], see also 
[Wu02, Ju05]. 

Lemma 2.6. Suppose that f,g€ C'“(H), then 

(2.2) l|A^(/ff)llL^ < c{\\A^f\Mg\\Lr2 + ||a*<7||l.i||/||l.2) 
where s > 0, 1 < r < Pi,P 2 ,qi,q 2 < oo and 1/r = 1/pi + l/p2 = l/o'i + I/q' 2. 

Lemma 2.7. Suppose that f,g & (7^(12), then 

(2.3) ||A^(/ • V5) - / • (A^V5)||l2 < C(||V/|Upi||A^<7||lp2 + ||A7||l« ||V<7||l.2 ). 

where s > 0, 2 < pi,P 2 ,qi,q 2 < oo and 1/2 = 1/pi + l/p2 = l/o'i + 1/^2. 

Remark 2.1. We remark that the inequalities (2.2) and (2.3) in Lemmas 2.6 and 2.1 are 
also valid for those f and g belonging to certain Sobolev spaces which make the right-hand 
sides of (2.2) and (2.3) to be finite. 
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2.3. The main results. We now state the result about the existence of weak solution and 
global strong solutions of the 2D Boussinesq system (1.1). The proof involves the standard 
Galerkin approximation, some basic functional analysis theorems in [Tem84], the uniform 
estimates in Section 3, and continuity estimates in Section 4 . 


Theorem 2.1. Let 


Hi = G l2(D) : j 0dx = 0 


and 


H 2 = I u £ : V • u = 0, / uidx = / U2dx = 0 

I J D. J o, 


Suppose f £ H ^ and {9 q,uq) £ Hi x H2. Then, for any T > 0, there exists at least one 
weak solution of the 2D Boussinesq equations (1.1) in the following sense: 

^ uipdx— f ui^u-V(p)dx+i' [ (A“u)(A“(^)da: = f 9e2udx, = {ipi,(p 2 ) £ , 

J ^ J n J Q. 


dt Jn 

and 


^ f difdx — f 9{u ■'Vfj)dx + K f {A^9){A^'ip)dx = f ffjdx, Vip£C°°{Q). 
dt Jq Jq Jq j 


Moreover, 9 £ L°°(0, T; Hi) n 1.2(0, T; H^^) and u £ L°°(0, T; H 2 ) n L2(0, T; H“). 

Furthermore, if we assume that si, S 2 satisfy (1.3) and (1.4), (6*0 , wq) £ x H^^ 
f £ n LP for p > 2, then for any T > 0, the Boussinesq system (1.1) has a unique 

strong solution {u, 9) satisfying 


{9,u)£ C([0,T],H^i)nC([0,r],H^=), 
{9t,Ut) £ L2(0,r;H"i-l^) nL2(0,T;H"2-“). 


The main goal here is to prove the existence of the global attractor for the Boussinesq 
system (1.1) and we have the following theorem. 

Theorem 2.2. Assume that v > 0, k > 0, si and S 2 satisfy (1.3) and (1.4), and f £ 
n LP for p > 2. Then the solution operator {S{t)}t>o of the 2D Boussinesq system: 
S{t){9o,UQ) = {9{t),u{t)) defines a semigroup in the space x H*^ for all t G M+. 
Moreover, the following statements are valid: 

(1) for any {9 o,uq) £ H*i x , t i->- S(t){9Q,UQ) is a continuous function from 1R+ 
into H®! X H®2; 

(2) for any fixed t > 0, Sft) is a continuous and compact map in H®i x H®^; 

(3) {S{t)}t>o possesses a global attractor A in the space x H®^. The global attractor 
A is compact and connected in x H®^ and is the maximal bounded attractor and 
the minimal invariant set in x H®^ in the sense of the set inclusion relation. 


2.4. The global attractor. In order to prove the main result Theorem 2.2, we are going 
to utilize the abstract result from [Tem88, Ghapter I] about semigroups and the existence 
of their global attractors. One can also refer to [Lad91], [BV92], [HK06] for additional 
development of the theory of global attractor. 

Theorem 2.3. Suppose that X is a metric space with metric d{-, ■) and the semigroup 
{5(t)}t>o is a family of operators from X into X itself such that: 


THE GLOBAL ATTRACTOR OF THE 2D BOUSSINESQ EQUATIONS 


7 


for each fixed t > 0, S{t) is continuous from X into itself; 
for some to > 0, S{to) is compact from X into itself; 

there exists a subset Bq of X which is hounded, and a subset U of X which is open, 
such that Bq CU C X, and Bq is absorbing in U for the semigroup, that is for any 
bounded subset B dU, there exists to = to{B) > 0 such that 

S{t)B C Bo, Vt > to(^)- 

(Bo is also called the absorbing set of the semigroup in U). 

Then A ;= uj{Bo), the u-limit set of Bq, is a compact attractor which attracts all the bounded 
sets of U, that is for any bounded set B cU, 

lim dist{S{t)B,A) = 0, 

t —^~ 1“00 

where dist{Bi, B 2 ) := miy^B 2 d{x,y) is the non-symmetric Hausdorff distance be¬ 

tween subsets of X. Furthermore, the set A is the maximal bounded attractor in U for the 
inclusion relation and the minimal invariant set in the sense of S{t)A = A,\/t> 0. 
Suppose in addition that X is a Banach space, U is convex and 

(iv) for any x € X, 1S{t)x from M_|_ to X is continuous. 

Then A = u}{Bq) is also connected. 

If U = X, A is the global attractor of the semigroup {S{t)}t>o in X. 

Remark 2.2. We will carry out the proof of Theorem 2.2 by checking all the items in 
Theorem 2.3. We first show the -uniform estimates in Section 3.4, which implies items 
(ii), {Hi), then we check item {iv) in Section 4-1 by proving that {S{t)}t>o is continuous 
from M_|_ to x , and finally we check item {i) in Section 4.2 by proving that {S'(i)}i>o 
is continuous in x . 


(0 

{ii) 

{Hi) 


3 . Uniform estimates 


In the following, we denote by C a positive constant, which is independent of time t and 
of the initial data uq and 6q. The constant C may vary from line to line. 


3.1. and L^-estimate for 9. 


Proposition 3.1 (Existence of absorbing ball in ifi and for 9). Under the assumptions 
of Theorem 2.2, there exists t^ = tidl^olliz) > 0 , such that 


(3.1) 

\\9{t)\\L2<C, ytPtl, 

and 

rt+^ 

(3.2) 

j IIA^^II^ads <C, yt> tl 


Proof. Here and throughout this article, we let Ai be the first eigenvalue of A. Then by the 
results from [Ju05, Section 5.1], we obtain 




ll/lliA ii/iij. 

K^Xf) Ic2Af’ 


(3.3) 


wmiih < 
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which immediately implies the uniform bound (3.1) for some > 0 large enough. Further¬ 
more, taking inner product of (1.1)3 with 6 in L^(D) and integrating in time give 

pt+l \\f\\2 

(3.4) ||«(4 + l)||J,+„ / ||A49||i,d»<||e(t)||i, + MM. 

Jt 

Therefore, the time average estimate (3.2) follows from (3.1). We also deduce the equation 

(5.4) in [Ju05] that for all 2 < p < oo, we have 


(3.5) 


\\^it)\\Lp < e 




P\\f\\Lp\ P\\f\\LP 

KXf ) ^ KXf ■ 


The above inequality gives the uniform estimate and absorbing ball in for 9 whenever 
9o € LP(H) for all p G [2, cx)). □ 


3.2. L^-estimate for u. 


Proposition 3.2 (Existence of absorbing ball in for u). Under the assumptions of 
Theorem 2.2, there exists = t2(ll^o||L2 ) ||'^o||l2) > 0, such that 


(3.6) 
and 

(3.7) 


\u{t)\\L2<C, yt>q, 


pt+l 

\\A^u\\l2ds <C, yt> t*2. 


Proof. Taking the inner product of the equation (l.l)i with u in L^, since u is divergent 
free and {u ■ 'Vu,u) = 0, we have 

+ + = (^ 62 ,n) < ^\\A-^e\\l2 + ^\\A°u\\l2. 


Then, 

(3.8) 

Since u and 6 have mean zero, by Poincare’s inequality, we have 
(3-9) < 7T2 tII^IIl2- 


^\\u\\l2 + n\\A^u\\l2 < ^\\A °e\\l2. 


Now, integrating in time and using (3.3), we obtain in the case that ^ 
(3.10) 


\u 


L2 




< ^\\uo\\l2 + 


uXf 


g—_g—KA^^t 




II, ||2 Il/lli2\ , WfWh 


2 

L2 


and in the case that z^A^“ = ztA^^, 


(3.11) \\u{t)\\l 2 <e 


-uXpt 


2 

L2 




2 

L2 


Af; KAf(z.Af)2’ 
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which shows the uniform estimate (3.6) for some > 0 large enough. Furthermore, inte¬ 
grating (3.8) in time gives 

rt+l 1 

(3.12) \\u{t + l)\\l.+u ||A“rr||i.ds < \\u{t)\\l, + 

Therefore, the time average estimate (3.7) of l|A“ti|p follows from (3.1) and (3.6). □ 

3.3. L^-estimate for the vorticity w. In order to proceed to find a uniform H^-estimate 
for (9, u) in the next subsection, we also need a uniform L^-estimate for the velocity u. 
However, we could not simply multiply \u\p~^u to show the L^-estimate because of the 
troublesome pressure term. Here, we aim to show the uniform L7^-estimate for the velocity 
u by proving the uniform L^-estimate for the vorticity w, which satisfies 


(3.13) 


dtuj + u ■ Vlj + i/(—A)“a; = di9, 

u = 


where V-*“ = (—92,9i)- 

Proposition 3.3 (Existence of absorbing ball in for w). Under the assumptions of 
Theorem 2.2, there exists tg = tgdlwollL^) II^oIIl^) > 0, such that 


(3.14) 

ll‘^(^)lll,T W'^Wh^ < C, 

\/t > tl 

and 

/■*+! 


(3.15) 

J IIA^wll^^ds < C, 

yt > t*2. 


Proof. Taking the inner product of the equation (3.13)i with ui in we have 
^ ^ + n\\A^uj\\l2 = {di9,uj) < ^\\A^-°e\\l2 + ^\\A^uj\\l2. 


2 dt 


Hence, 

(3.16) 

Since uj has mean zero, Poincare’s inequality implies 

(3.17) 






By the assumption (1.2), we have 1 — a < /3 and using the interpolation inequality, we 
deduce from Proposition 3.1 that 

\\A^-^9\\l2ds <C, Vt > tl. 


I 


Hence, applying Lemma 2.4 with y = ||w||^ 2 , 9 = and A = we obtain 

the uniform estimate (3.14) for w and also for u by the Biot-Savart law. Furthermore, 
integrating (3.16) in time yields 


rt+^ 1 /-t+i 

(3.18) Mt + I)\\l2+n \\A-co\\l2ds < Mt)\\l2 + - / IIA^-^^Hi^ds. 

Jt ^ Jt 

Therefore, the time average estimate (3.15) of ||A“w|p follows. 


□ 
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As an immediate consequence of Proposition 3.3 and the Sobolev embedding theorem, 
we have a uniform L^-estimate for u and a time average estimate of ||A^'’""ti|p, that is for 
all 1 < p < oo, 

(3.19) II'^IIlp < C{p), j ||A^+“u||^ds < C, Vt > tg, 

where the constant C{p) >0 only depends on p, but is independent of the time t and the 
initial data uq and Oq. 

3.4. H^-estimate of 9 and u. 


Proposition 3.4 (Existence of absorbing ball in for {9,u)). Under the assumptions of 
Theorem 2.2, there exists t\ = t 4 (||^ollL 2 ) II'“oIIl2) > 0; such that 


(3.20) 

and 


l|0(t)lkn <G, 


|w||//S2 < C, Vt > tl, 


rt+1 rt+l 

(3.21) \\A^^+f^e\\l,ds<C, ||A 


S2+Oi.,\\2 


u\\ T2ds < C, 


\/t > t 


^ ('4- 


Proof. Taking the inner product of the equation (1.1)3 with in we have 


1 d 


(3.22) 


--^e\\l2 + {u ■ ve,A^^^9) + K||A*i+^0||i, = (/, a2^i0) 


2dt 




Since u is divergence free, then n • V0 = V • {u9) and according to Lemma 2.6, we have 


(3.23) 

\{u ■ V0, A2"i 0)| = |(V • (u0), A2 *i0)| 

< ||A*i+^i0||i2||A^i+i-^i(«0)llL2 

for some /3i > 0, pe, qg > 2, which are determined later, such that /3i < /3 and l/pg + l/qg = 
1/2. In order to determine pg, qg and /3i, we first choose ri > 0 such that 


(3.24) ri 


si, 2max{l — a, 1 — /3} < si < 1, 

any number in (2max{l — a, 1 — /3}, 1), si > 1. 


Then, the Sobolev embedding theorem implies that 


00 e C C L?®, 


where qg is chosen such that 

1 1 — ri 

qe 2 


1 

2’ 


< min{a, /3} 
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(3.25) 


Therefore, by the uniform RP-estimates (3.5) and (3.19) for {u,6), we infer from (3.23) that 
for t > max{t);, 

= : h + h. 

Now, we set si + 2 — 2/p0 — /3i = si + (3i, that is /3i = 1/2 + I/qq < min{a,/3}, then the 

Sobolev embedding theorem implies that pe ^^{Q) ^ and the term 

Ii can be estimated as 

h < C||A"i+^i0||i2||A"i+^i'u||i2 

(3.26) ^ ^ 

where we used the interpolation inequality for the tuples (si, si + /3i, si + /3) and (si, si + 
/3i, Si + a). Since 0 < si < S 2 , applying Poincare’s inequality, we find 

(3.27) ||A"i+"-u||i2 < C||A"2+“u||^2, \\A^^u\\l2 <C\\A^^u\\l2. 

Hence, using Young’s inequality, we find 

h < C||A*i+/^0|||||A^i0||^;^||A^2+a^||f ||A^2«||^;^ 

< + §l|A“ 9 |li. + jl|A-+”«|li. + ^l|A««||i., 

where a = f3i / {/3 — Pi) and b = Pi/{a — Pi). For the term I 2 as for R, we deduce from the 
same Sobolev embedding and interpolation inequality that 


(3.28) 


(3.29) 


2/^1 n 2/^1 

h < C||A"i+^i 0||22 < C||A"1+^0||J ||A*i0||j ° 


Il2 


C 


<Y^I|A“'+'’9|li, + ^||A"9|||,. 


Finally, we arrive at a differential inequality for 6: 

l^iiA*‘9iiia + <jiiA*‘+'’«iii2 + jiiA«+“«iii, + T||a*.«i|2, 

* * +fl|A"-'’/lli. + ^l|A»u|li,. 

Taking the inner product of the equation of (1.1) 1 with hP^'^u in Lp, we have 
(3.31) ~^\\A^^u\\l, + {u -Vu.A^^^u) + u\\A^^+^u\\% = {ee 2 ,A^^^u). 

The right-hand side of (3.31) is estimated by 

\{ee2,A^^^u)\ = |(A"i+^0e2,A2*2-si-/3.j^^| < \\A^^+Pe\\L2\\A^^^-‘^^-f^u\\L2, 

which by using the interpolation inequality for the tuple ( 0 , 2s2 — si — P,S2 + a) because of 
the assumption 0<S2 — si<q;-|-/ 3 , is bounded by 


||A^i+/'0||^2||A 


S2+a 


^ 3i+a+l3-S2 


It 


1,2 


S2+“ 


I It 


si+a+^ —«2 

S2+“ 

L2 : 
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and finally applying Young’s inequality, and since ||■u|| 2,2 is uniformly bounded by Proposi¬ 
tion 3.2, we obtain that for t > 

(3.32) K«e2,^||A«+“u||i, + ^U-r, 

where a' = (si + a + /3 — S 2 )/(s 2 + a). We now estimate the nonlinear term in (3.31). 
Similar to (3.23) and using the L^-estimate (3.19) of u and Lemma 2.6, we first have that 
for t > max{t 2 ,t 3 }, 

I {U • Vu , I < C|| A"2+"lii||^2 IIyV"2+l-ai ^ 11^2 

(3.33) < C||A"2+“iii||^2 ||A*2+i-“iit||iP„ \\u\\l<,^ 

< C||A"2+“iw||i2||A'*2+i-“iu||iP„, 


for some ai > 0, Pu, Qu > 2, which are determined later, such that ai < a and 1/pu + 
l/O'tt = 1/2. Now, if we set S 2 + 2 — 2/pu — ai = S 2 + «i) so that ai = 1/2 + 1/qu = 
1 — 1/pu and choose large enough such that ai < a, then by the Sobolev embedding 

theorem U* °^^(^) the interpolation inequality for the tuple 

(s 2 j S 2 + 'S 2 + «) and Young’s inequality, we arrive at 


(3.34) 


20.1 

\{u ■ Vu,A^^^u)\ < C||A" 2 +«iii ||22 < C||A" 2 +“ii||^ 

<r !lll A^2+a„,||2_ , ^11 A^2„,||2_ 


||A 


S2 


u\ 


1 , 2 -- 

L2 


where b' = ai/{a — ai). Hence, we infer from (3.31) that 
(3.35) 

1 d ,,o ,, . ,,9 K,, . ,o V 


2dr' 


< jllA* 


2 

L2 


c 


+ -||W^+“r.||i. + ^||A 


u 


|2 _L 

Il2 + 


c 


! 


-1 


Summing the differential inequalities (3.30) and (3.35) together, we obtain that for t > 


(3.36) 




< 


C (4;l|A-9|li, + (4 + L)l|A««||t,) + ^||A--'> 


/II 


2 4 - 
L2 + 


c 


K'^ Vo.' 


T-1 


Here, starting with si = > 2max{l — a, 1 — /3}, < /? and S 2 = = 1, which 

implies that < min{2a — 1, 2/? — 1} < a + /3, then by Poincare’s inequality and 

equations (3.2), (3.19), for t > t* := max{t^, ^ 3 }) we have 


(3.37) 


/■i+l (1) ^t+l 

/ ||AN e\\l,ds < C, / ||A 

Jt Jt 


.( 1 ) ,,9 

’2 ull^ads < C. 


Hence, applying the uniform Gronwall lemma to (3.36) and using (3.37), we obtain that for 

t > r + 1, 

||A^“0(t)||L2<C, ||W2^’r.(t)||i. <C, 


(3.38) 
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Furthermore, integrating (3.36) and using (3.38), we find that for t > t* + I, 

(3.39) / ||A"i +^e\\l2ds <C, / ||A"2 +^u\\l2ds < C. 

Jt Jt 

Then, we iterate with si = + (3 and let S 2 = s® be any number that satisfies 

^2^^ — < a + f3 and + a. Hence, we have from (3.39) that for t > t* + 1, 

(2) ft+i 12) 

(3.40) / ||AH e\\l2ds<C, / ||A"2 u\\l2ds<C. 

Jt Jt 


Applying the uniform Gronwall lemma to (3.36) again, this time using (3.40), we obtain for 
t>t* + 2, 

(3.41) \\A^'"'eit)h2<C, \\A^'^'umL^<C, 


and 

(3.42) 


l'^'\\A^i'+Je\\l2ds<C, 



( 2 ) 

IIA'*^ +“it||i2ds < C. 


Therefore, with a bootstrapping argument, for any given real numbers si,S2 that satisfy 
Si > 2max{l — a, 1 — /3}, S2 > 1 and S2 — si < a + /3, (3.20) and (3.21) are proved. In 
addition, for fixed T > 0, 

(3.43) r ||A"i+^0|||2ds < oo, r l|A"2+"-u|||2ds < oo. 

Jo Jo 


By (3.41), given (6q,uq) € x then for some t > 0 large enough, the solution 
{6{t), u{t)) of system (1.1) belongs to the space x H*"* for some S3 = si + /3 and some S4 
such that 0<S4 —S3<a + /3 and S2 < S4 < S2 + a. By the Sobolev compactness embedding 
theorem in [Tem88], the inclusion map x H*"* i-;- x is compact. Thus, for any si 
and S2 satisfy (1.3) and (1.4), the solution operator S{t) defined by S{t){6o,uo) = {6{t),u{t)) 
is a compact operator in the space x for some t > 0 large enough. □ 


4. Continuity 

4.1. Continuity with respect to t. In this section, we check that the solution operators 
{S'(t),Vt > 0} are continuous in the space x with respect to t. In order to achieve 
this, we recall the following lemma which is a particular case of a general interpolation 
theorem in [LM72]. A proof can be found in [Tem84]. 

Lemma 4.1. Let V, H, V' he three Hilbert spaces such that 

V C H = H' cV 

where H' is the dual space of H and V is the dual space ofV. 

If a function u belongs to L^(0, T;U) and its derivative u' belongs to LP‘{d,T;V'), then u is 
almost everywhere equal to a function continuous from [0, T] into H. 

Proposition 4.1. Under the assumptions of Theorem 2.2, the solutions of Boussinesq 
system (1.1) satisfy {9,u) & C{[0,T]; x C([0, T]; H®^). 
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Proof. For any fixed T > 0, the uniform estimates proved in Proposition 3.4 ensure that 

{9,u) G X L2(0,r;H"2+"), 

where si and S 2 satisfy (1.3) and (1.4). Hence 

(A"i0,A®2-a) G X L‘^{0,T; H^). 

Our goal is to show that {0,u) G C{[0,T]; x C{[0,T]-, H^^), that is {A^^9,A^^u) G 
C([0,T];L^) X C([0,T];L^), and by Lemma 4.1, it suffices to show that (A^^^*, g 

L^([0, T]; H“^) X L^([0, T]; iL“"). For any {ip,(p) G x (iL“)^, we derive from (1.1) that 


(4.1) 


' {A^^9t,fj) + {A^^{u ■V9),f;) + K{A^^+^f^9,fj) = (A^i/,V'), 
{A'^^ut, (f) + (A^2(ii • Vu), (f) + (p) = (A®26»e2, ^p). 


Then, by the Cauchy-Schwarz inequality, 

(4.2) 

fKA^i0i,V')| <||A^i-/'(tx.V0)|U2||V^||^,+K||A^i+/^0||i2||V’|b/3 + l|A^^-^/llL2||V>||j^/3, 
||(A^2^t,(^)| < \\A^^-°‘{u ■ Vu)\\l2\\(p\\h<- + iy\\A^^+°‘u\\L2\\(p\\H<^ + \\A^^-°‘9 \\l2\\(p\\h<-, 
which implies that, 

f WA^^OtWu-^ < ||A^i-/^(ti ■ V0)||i2 + K\\A^^+f^9U2 + ||A*i-/^/||i2, 


(4.3) 


WA^^utWn-c^ < ||A^ 


[u 


Vu)\\l2+i^\\A^^+^u\\l2 + \\A^^-<^9\\l2. 


We now estimate the nonlinear term u-'S/9 and ii ■ Vit in the right-hand side of (4.3). Since 
u is divergence free, we have 

(4.4) ||Wl-^(^l • V9)\\l2 = ||Wi-/^V • {u9)\\l 2 < \\A^+^^-^{u9)h2. 

We now choose r > 0 such that 

si, 2max{l — a, 1 — /3} < si < 1, 

any number in (2max{l — a, 1 — /?}, 1), si > 1. 


(4.5) 

Then, 

where 


r = 


9o G CH^ C LP\ 


1 


1-^ r 1 1. 1 

— <mm{a-2.<9-2}<5 


Uq G cH^ C LP\ 
2 


P2 = 


Pi 2 ' 2 2 2 ^ 2/3—1 

Hence, by (3.5) and (3.19), 9 G L°°([0,-|-oo); L^i) and u G L°°([0,-boo); L^^). Therefore 

(4.6) \\A^+^^-^iu9)\\L2 < +C||tx||LP.||Ai+^i-/'0||i,., 


where 1/pi -b 1/g'i = 1/2 and l/p 2 + 1/^2 = 1/2. Now, we choose ql such that 

1 S2 “b (U ~ (1 “b Si — /3) 1 

7i ^ 2 = 2 - 

Since 1/qi = rjl and S 2 > si, we have 

4<l(2-a-^)<^ - 2max{l - a, 1 - ^} < = —. 

2 2 2 qi 
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Hence, qi < qX- Thus, by Poincare’s inequality and the Sobolev embedding theorem, 

(4.7) ||A^+"i-^-u||l.i < < G||A"=+“it||i2. 

We also deduce from l/q 2 = 1/2 — l/p2 = 1-/3 that 

1 Si + /I — (1 + Si — /3) _ 1 

^ 2 “ 2 ’ 

and hence by the Sobolev embedding theorem, 

(4.8) < C||A^i+^0||i2. 

Hence, by (4.2), (4.4), (4.6),(4.7), (4.8), we have 

<C||0|U.i||A*^+“ti||^2+C(||tr||iP.+K)||Wi+^0||i2+C||Wi-^/||i2, 


which by utilizing the estimates (3.5), (3.19) and (3.43), shows 

rT 


f 


IIA^i^tll^.^ds < oo. 


Similarly for the term u Vu, we similarly have 

||A*2-“(wVtr)||i2 = ■ {u(^u)\\l2 < ||A^+^2-“(tr ® w)||i2 


(4.9) 


< C'||u||i,p3||A^+^2 °‘u\\li3, 


where l/ps + l/gs = 1/2 and q^ = 1/(1 — a). By (3.19), we know that u € L°°(0, +oo; 
for p 3 = 2!{2a — 1). Since 

1 ^ S2 + ct — (1 + S2 ~ cr) 1 
® ^ 2 ^ 2’ 
then the Sobolev embedding theorem shows that 

(4.10) ||A1+"=-"w||l.3 < C||A"2+“-u||^2. 

Therefore, by (4.3), (4.9), (4.10) and the assumption that 0 < S2 — si < a + /3, we find 

(4.11) WA^^utWn-c. < Ci\\u\\LP3 + z.)||A*2+“ri||i2 + C||A^i+^0||i2, 
which by combining the estimates (3.19) and (3.43), yields 

fT 

/ WA'^^utWjj-cds < oo. 

Jo 

We thus completed the proof of Proposition 4.1. □ 


4.2. Continuity for fixed t > 0. We simultaneously prove uniqueness and continuity of 
S{t) from x to itself for any hxed t > 0. 

Proposition 4.2. Under the assumptions of Theorem 2.2, the solution of Boussinesq sys¬ 
tem (1.1) is unique and the solution operator S{t) : x 77*^ i—)■ 77*^ x 77^^ is continuous 

for any fixed t > 0. 
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Proof. Let (ui,0i,vri), {u2,92,'^2) satisfy (1.1) with two initial data {u^, 0?) re¬ 
spectively. Then Q = ui — U 2 , rj = 9i — 02, 'k = — 7^2 satisfy 

f dtC + ui - VC + C ■ Vu 2 + z^(-A)"C = -Vvr -|- rje 2 , 

\dtr] + ui -Vr] + C ■ V02 + K{-A)^r] = 0. 


(4.12) 


Taking the inner product of (4.12) with (A^^^^^^ A^^^r/) in and using V ■ A^^i^ = 0, which 
comes from the fact that A* and V commute with each other and is divergence free, we 
obtain that 
(4.13) 

I + ^^l|A^^+“Clli2 = {ve2,A^^K) - {ui ■ VC, A^^^o - (C ■ Vti2, A^^^C), 

Usl|A^^??lli2 + K||A^i+^r?|||2 = -(ui • Vr/,A 2 ^ir/) - (C • V6I2, A^^ir/). 

Similar to (3.32), using Poincare’s inequality, the interpolation inequality and Young’s in¬ 
equality, we have 

I(r?e2, A^^^OI < ?I|A*^+Mli. + ^||A*^+“Clli2 + ^ 


(4.14) 


-I 

6 ' 

K , 


, 


<^l|A"+Mi!« +jiiA«+”ciii.+ 


c 


-1 


2 

L2 


1/ a' 


-llA^^ciii2, 


To deal with the term {ui ■ we observe that {ui ■ VC)^ C) = 0; so that 

(4.15) {ui ■ VC, A^^^C) = (A*^(«i • VC), A^^C) = {A^^ui ■ VQ - • V(A*2C), A^^C)- 

Noticing that V and A commute, we find 

|(A^=(ni • VC) - u, • V(A^X), A^2^)| = \{A^^{ui • VC) - u, • (A^^VC), A^^OI 

< CWA^H^^i ■ VC) - u, • (A^^VC)k2||A*^CIlL2. 

Applying Lemma 2.7 for 
(4.17) 


, 11111 

Pi,P 2 ,qi,q 2 > 2 and- ^ -=- ^ -=-, 

Pi P 2 qi q 2 2 


we have 
(4.18) 
Let 


||A^2(rii • VC) - ui • (A^2VC)||l2 < CiWVuiWLPi W^C\\Lr 2 + ||A^2^i||l,i ||VC||l.2; 

< CdlAuilliPi ||A®2 C||lp 2 -I- IIA^^uilli,?! ||AC||l'J2) 


Pi = 


P2 = 


qi = 


q2 = 


a ^ 1 — a 1 — a 

Then, using the Sobolev embedding inequalities, we obtain 

IIAuilliPi < C||A2 -“ui|Il 2 < C||A^2+“txi||i2, 


a 


||a^^CIIlp2 <c||a*2+“CIIl2, 

||AC||L.2<C||A2-“C||i2<C||A^^+“CllL2. 
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Then, 

(4.19) ||A^=(t.i • VC) - • (A^=VC)||l2 < C||A^=+“t.i||^2||A*2+“C||L- 

Hence, using (4.15)-(4.16) and Young’s inequality, we have 

1(1.1 • VC, A^^^OI < C||A«2+“i.i||i2||A^^+“C||L2||A^^CIlL2 


(4.20) 


<^||A^^+“t.i||i.||A^^CIli2 + ^l|A*^+“Clli2. 


By the Cauchy-Schwarz inequality, we have 

KC • Vi.2, A^^^C)! = l(A*^-“(C • ViX2), A*2+“C)I 


(4.21) 


<^||A*-“(C-VtX2)||i2 + ^||A^^+“Clli2. 


Again we apply Lemma 2.6 for pi, p 2 , qi, q 2 satisfying (4.17) to obtain 
(4.22) ||A^^-“(C • Vt.2)||L2 < C(||A*2-“C||l.i ||Ai.2||l.2 + ||C||l^i ||A*2-“+1ix2||l.2). 

Let 


1 


Pi = 


P2 = 


qi = 


q2 = 


1 — a a 2a — 1 

Then, using the Sobolev embedding inequalities, we have 

||a*^-“CIIlpi <c||a*^CIIl2, 


||At.2||LP2 < C||A2-“t.2||L2 < C||A^^+“iX2||l2, 

IICIIl.1 <C||A2-2“C|L2<C||A^^CIIl^, 

||A*^-“+'i.2||l.2 <C||A^^+“i.2||l2. 

||A^^-“(C-Vi.2)|li2<C||A*^CIli2||A^^+“ 


1 — a 


and 

Then, we have 

(4.23) 

and hence, 


"“211^2, 


(4.24) 


KC • Vi.2, A^^^C)! < ^||A^^CIli2||A^^+“t.2||i2 + pA«^+“Clli2. 


Next, we estimate the term {ui ■ Vrj, Since (iti • V{A^^rj), A^^r]) = 0, we have 

(4.25) (wi • Vp, A^'ip) = (A"i(i.i • V 77 ), A^ip) = (A"i(i.i • Vr?) - ui ■ V{A^^r]),A^^rj). 
Again, since V and A commute, we have 

I (A"i (wi • V 77 ) - ui ■ V(A*i rj ), A"ir?) I = I (A*i {ui ■ Vr?) - ui ■ (A"i Vr?), A*ir?) | 

< C-IIA^H^i • Vt?) - 1.1 • iA^^Vrj)h2\\A^^rj\\L2. 
Applying Lemma 2.7 for pi, p 2 , <?i, q 2 satisfying (4.17), we have 

\\A^^{ui • Vr?) - 1.1 • (A^iVr?)||i 2 < C(||Vi.i|Upi \\A^^v\\lp 2 + ||A*1i.i||l,i ||Vr/l|L«) 

< CdlAitillLPi ||A®ip||i,P2 + 1|A^1 i.i||lpi ||Ar?||LP2). 


(4.27) 
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We set Pi = 2//3 and p 2 = 2/(1 — j3). Then, using the Sobolev embedding inequalities, we 
obtain 

IIAwilliPi <C\\K^-Pui\\l2<C\\K^^+^ui\\l2, 

and 

Let 7 = min{a,/3}, qi = 2/(1 — 7) and q 2 = 2 / 7 . Since si > 2 max{l — a, 1 — /3}, then 
Si > 2(1 — 7). Thus, we have 2 — 7 < si + 7, and 

||A^itxi|U.i < C||A^i+^txi||i. < C||A*i+“txi||^2 < C||A^^+“tii||i2, 


||Ar?|U.2 < C||A2-7r?||i. < C||A*i+^r?||i. < C||A*i+M|i2. 

Hence, 

(4.28) ||Wn«i • Vr?) - til • (A'*^Vr?)||i 2 < 

Therefore, using (4.25)-(4.26) and Young’s inequality, we have 

|(ui-VC.A2*‘,)|<C||A«+"u,||i.||A-+7lli>l|A*‘.,||y 

< ^l|A«+“t.,||i,||A"^||i, + b|A«+7[||,. 

K D 

For the term {(^ ■ V02)A^*ir/), we consider two cases: si > /3 and si < /3. Suppose that 
•si > A, so that 

(4.30) |(C-V02,A2*ir?)| = |(Wi-^(C-V 02),A^^+^)| < ^||A^i-/3(C.V02)|li2 + ^||Wi+^||i2. 

tx O 

Applying Lemma 2.6 for pi, p 2 ) <?i, <12 satisfying (4.17), we have 

(4.31) ||a^i-/^(c-V02)||l2 <C’(||a^i-/^CIIlpi||A02||lp2 + ||CIIl^i||A"i-^+102||l^2). 

Let 


1 


Pi = 


1-/3’ /3’ 2,5-1 

Then, using the Sobolev embedding inequalities, we obtain 

||a^^-^CIIl^i < cwa^KWl^ < c||w=CIIl2, 


qi = 


q2 = 


1-/3 


and 


||A 021 |lp 2 < C\\A‘^-^ 92 \\l 2 < C||A*i+/^ 02 ||ll 

IICIIl.1 <c\\a^-^^c\\l^<c\\a^x\\l 2 , 

||A^i-^+102||l<^2 < C||A^i+/^02||l2. 

Then, we have 

(4.32) ||Wi-^(C-V02)|li2 <C||A*^CIli2||Wi+/'02|li2. 

Hence, using Young’s inequality, we have 


(4.33) 


|(C • V02, A2*i7)l < -||W^CIli2||W^+^02|li2 + ^||Wi+/^7lli2. 

Kj o 
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Otherwise, si < /3 so that 2si < si + /3. Using the Cauchy-Schwarz inequality, Poincare’s 
inequality and Young’s inequality, we have 

|(C-V02,A2 *ir?)| < ||C-V02 ||L2||A2^ir?||^2 

<C'||C-V02||L2||A^^+M|i2 


(4.34) 


C, 


<-llc-v02iii. + wi|A^^+Mli2. 

K O 

For any pi,p 2 > 2 such that l/pi + l/p 2 = 1/2, we have 

||C-V02||l2 < ||CIIlpi||V02||l.2 <C||C||l.i||A02||lp2. 


Let Pi = 2/(1 —/?) andp 2 = ‘^IP- Since/3 < 1 < S 2 , using the Sobolev embedding inequality, 
we have 


(4.35) 
and 

(4.36) 
Hence, 

(4.37) 


IKIIlpi <c||a/^CIIl2<c||a^^CIIl2, 

||A02||lp2 < c\\A^-^e2h2 < c||a*i+/'02||l2. 


KC • V 02 , < ^||A*^CIli 2 ||A*^+^ 02 |li 2 + 


Therefore, using (4.14), (4.20), (4.24), (4.29), (4.33) and (4.37), we have 
(4.38) 


i4||A“>cili2 + i'liA»+“ciii> <^I|A«+Mii. + ^iiA«+“ciii2 


+ C(||A«+“,.,lli, + ||A«+“,.2||i,)||A‘»<||i„ 


and 

(4.39) 


14||A*>,[|i,+K||A"+M|i, < |||A«+'’()|||, 


+c’iiA«+“.*iiri2iiA«,iri2+ciiA«+'’e2iii2iiA«<iri2. 


Summing (4.38) and (4.39) gives 

4{||A*xiii. + liA*‘i(iiio + ■'I|a*>+“cii! 2 + »=I|A*‘+Mii, 

< C(||A«+”ui|| 2, + ||A‘>+“u2lli2 + l|A”+'’«2lli2)l|A»C||j2 + 

< C(||A«+“ui||= 2 + ||A«+“u2|li2 + IIA-'+^OjIII^XIIA'^CIIJ, + ||A'‘„|||,). 

and using the Gronwall inequality in Lemma 2 . 2 , we obtain that 

ilA»c(*)lii> + liA"x(*)lii. 

<C(||A«C(0)|li2 + ||A",(0)||i2)exp ||A«+“ti,(s)l||, + ||A*-+“u 2 {s)||i 2 + ||A«+'’92(o)lli.ds 

Notice that from (3.43) that 


f ||A® 2 +“it^(s)|| 22 ds < 00 ,/ ||A^’ 2 +"u 2 (s)||i 2 ds < oo, and / ||A*i+^ 6 ' 2 (s)||i 2 ds 

Jo Jo Jo 


< oo. 
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By the Riesz lemma, since ||A®ir/(0)||i2 and ||A^2 ^(q)||^2 go to zero, then \\A^^r]{t)\\i 2 and 
||A®2^(t)||^2 converge to zero for almost every t. Since ||A®i7/112,2 and ||A*2^||^2 are continuous 
in t as we have proved in Section 4.1, then (||A®ir/(0)||2,2, ||A®2(^(0)||2^2) converges to zero for 
every t. □ 
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